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An analytic method is given for deriving the part of the retarded Green's function 

\ v(x,x') that contributes to the tail term in the radiation reaction force felt by a 

O ' 

{Sj ' particle coupled to a massless minimally coupled scalar field. The method gives an 



expansion of v(x,x') for small separations of the points x, x' valid for an arbitrary 
static spherically symmetric spacetime. It is obtained by using a WKB approxima- 
tion for the Euclidean Green's function for the massless minimally coupled scalar 

t^j- ■ field and is equivalent to the DeWitt-Schwinger expansion for v(x,x ). The first few 

CO ' 

. terms in this expansion are displayed here for the case of Schwarzschild spacetime. 
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This latest version has two new appendices A and B containing respectively the two Errata 
in Ref. 153] and 15$ to the version published in Physical Review D f5£ j. In this version we 



^ ■ have made the corrections indicated therein. 

H ; 

The past several years have seen a resurgence of interest in the gravitational radiation 
reaction problem in part because space based gravitational wave detectors such as the Lasar 
Interferometer Space Antenna (LISA) are expected to be able to detect the gravitational 
radiation emitted when a compact object such as a stellar mass black hole or neutron star 



spirals into a supermassive black hole 



|2j. In such a situation the change in the spacetime 
geometry of the supermassive black hole due to the presence of the compact object and 
the gravitational radiation it emits is negligible (except near the compact object) so it is 
sufficient to compute the trajectory of the compact object in the background geometry of 
the supermassive black hole [2]. The radiation reaction under these simplifying conditions is 
known as the "self-force" . 

The radiation reaction problem for a point charge radiating electromagnetic waves in 
a curved space background was first investigated by DeWitt and Brehme [3] and later by 



Gravitational radiation reaction for a moving particle was considered by Mino 



Hobbs ^ 

et al ja] and Quinn and Wald [6]. The self-force on a particle interacting with a massless 
minimally coupled scalar field was considered by Quinn [7|. For a review see Ref. js]]. The 
scalar field case is of interest primarily because it contains many of the same features as 
the gravitational and electromagnetic cases while being much easier to work with. As found 
originally by DeWitt and Brehme ^3j for the electromagnetic case, it is always possible to 
break the radiation reaction force into two parts, one that is local and which reduces to the 
standard Abraham-Lorentz-Dirac expression in the flat space limit, the other part which is 
nonlocal and consists of an integral of the retarded Green's function over the past trajectory 
of the particle 5|, |6|, |7J • The nonlocal part is often called the "tail" term. 
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In most schemes that do not involve the weak field limit it must be computed numerically 
by expanding either the retarded Green's function, or the field itself, in terms of spherical 
harmonics in Schwarzschild or the spin-weighted spheroidal harmonics in Kerr [k], 
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29] . However, a drawback of these methods is that it has 



not been possible to separate out the local and nonlocal parts. The result is that there are 
divergences which must be regularized. The usual way to treat these divergences is through 
subtractions that occur at the level of the modes, resulting in a finite mode sum. This is 
similar to adiabatic regularization 



field theory in curved spacetime 




34| introduced in the context of quantum 



rzl 



361 ] . Zeta function regularization has also been used 
There is another way to obtain the tail term, at least in principle, by computing the part 
of the retarded Green's function that contributes to it. In this paper we propose to use the 
Hadamard expansion of the retarded Green's function for this purpose. It is known that 
this expansion is valid only when the points are close together. The self force on the object 
at a particular spacetime point has contributions from all points over its past trajectory. 
In nonlocal processes it is not unreasonable to start by taking into account contributions 
from points closest to the object as they usually give a greater weight than those farther 
away. A quasilocal expansion such as the one we are suggesting is the logical way to start 
an analytic approximation. Even if it turns out not to capture the dominant contribution, 
it is still worthwhile to investigate its range of validity for the specific task at hand. If any 
such analytic approximation produces even marginally reasonable results, it can provide a 
relatively quick way of estimating the self-force for a given trajectory. If the results are 



3 



accurate enough then it can be used in place of brute force numerical integrations, thus 
giving an immense economy of effort. 



For a massless minima 
Hadamard form 
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ly coupled scalar field, the retarded Green's function takes the 
3^ 1 



G ret (x,x') = 6(x,x') {^^l8[a{x,x')\ - ^le[-a(x,x')}\ (1.1) 

^ 47T 57T J 

Here 9(x, x') is defined to be zero outside of the past light cone and one inside of it. The 
quantity a(x,x') is equal to one-half the square of the proper distance between the points 
x and x' along the geodesic connecting them. The function v(x, x') contributes to the tail 
part of the self-force. It obeys the equation 



D x v(x,x') = . 



1.2) 



Note that v(x,x') is finite when the points come together. 

For small separations of the points the functions u and v can be computed by using the 



DeWitt-Schwinger expansion 
the result for v is 
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42]. To second order in derivatives of the metric 



v (x, x') 
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Christensen's expansions 41[ can be used to compute v(x, x') to fourth order. Phillips and 
Hu 39j have computed it to sixth order for a conformally coupled massless scalar field. Their 
result is useful for the study of a massless minimally coupled field in Schwarzschild spacetime 
because in any spacetime with zero scalar curvature the Green's functions for the massless 
scalar field do not depend upon the coupling to the scalar curvature. However, taking the 
general expressions for the terms in an arbitrary spacetime from Ref. 39] and evaluating 
them for the Schwarzschild geometry is still nontrivial. 

The results of these calculations show that for Schwarzschild spacetime the sixth order 
contribution to v(x, x') from the DeWitt-Schwinger expansion gives the first nonvanishing 
term, and that is the only term which the results from Ref. [39J] can provide. To compute 
the self-force more terms are needed from an expansion for v(x, x'). Also, given the fact 



1 Our definitions for the Hadamard expansion are equivalent to those of Ref. [37|, [38|, |39| . The conventions 
of Ref. 0] differ from ours. To obtain the convention for G rc t in Ref. 0] let G re t — * G re t/(47r). Substituting 
this into Eq. (|A1[) and comparing with Eq. (9) of Ref. 0] gives u(x, x') = U(x,x'), and v(x, x') — 2V{x,x'). 
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that no explicit computations of v(x, x') have been done for Schwarzschild spacetime us- 
ing the DeWitt-Schwinger expansion, it is of some general interest to obtain the first few 
nonvanishing terms in the series. 

In this paper we describe a method we have developed to compute an expansion for 
v(x, x') for a massless minimally coupled scalar field in a general static spherically symmetric 
spacetime. The method involves the use of a WKB expansion for the radial mode functions of 
the Euclidean Green's function for the scalar field. It gives the same results as the DeWitt- 
Schwinger expansion. Because of its reliance on a WKB expansion for the radial modes 
of the Euclidean Green's function, it should be straight-forward to generalize the method 
to the cases of the electromagnetic and gravitational fields in static spherically symmetric 
spacetimes. Since the scalar, electromagnetic, and gravitational wave equations in Kerr 
spacetime are separable 43, 44, 4^], it is likely that the method can be adapted to those 



cases as well. 

In Sec. [Til we show the relationship between the Euclidean Green's function and v(x,x') 
for the massless minimally coupled scalar field in a general static spherically symmetric 
spacetime. In Section [TTTJ our method of deriving an expansion for v(x, x') in powers of 
(x — x') is given and the first few terms of the expansion for Schwarzschild spacetime are 
displayed. We draw our conclusions in Section [TV] 

II. EUCLIDEAN GREEN'S FUNCTION AND THE TAIL TERM 

The metric for a static spherically symmetric spacetime can be written in the form 

ds 2 = -f{r)dt 2 + h{r)dr 2 + r 2 dQ 2 . (2.1) 

where f(r) and h(r) are arbitrary functions of the radial coordinate and dVt 2 is the metric of 
a 2-sphere. Changing the time variable t to t — it gives the Euclidean metric. (Quantities 
defined in the Euclidean space carry a subscript E). The Euclidean Green's function is a 
solution to the equation 

D x G E (x,x') = - 6 -^^- (2.2) 



where g E is the determinant of the Euclidean metric. If the analytic continuation 

(r-r') 2 -> -{t-t'f + ie (2.3) 



is used with e an infinitesimal positive real quantity then 46 1 



Ge(— it, x; —it', x') = iGp(x, x') (2.4) 



with Gp the Feynman Green's function. Since 



ImG F {x,x') = -^G {1 \x,x') 
where G^(x, x') is the Hadamard Green function, it is clear that 



(2.5) 



ReG E {-ir,x; -ir',x') = -G {1) {x,x') = .., 

2 87T^CT(X, X'J 107T" 1 



M (x,x') f ^^ log[|CT(X;X , )|] + ^^0 (26) 
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where in the second equality we have expressed G^(x, x') in the Hadamard form 37l.l38l.l39l]. 
Note that if the points are separated in the time direction, or if they are farther apart in 



the time direction than in other directions then 

1 
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a(x, x') = --f(r)(t - t'Y + 0{{x- x') s } 



(2.7) 



Thus v(x,x') is proportional to the coefficient of the log[(r — r') 2 ] part of Ge{x,x'). It is 
this fact that will enable us to derive an expansion for v by finding one for Ge using a WKB 
approximation. 

To find such an expansion for Ge first note that in a spacetime with metric (12. ip an exact 



expression for Ge is 



47] 



G K (x,x') 



poo 00 

— / (Iujcos[uj(t - t')} ^(21+ l)Pi(cos7)C^i£w(r<) qu(r>) 

71 ^° E=0 



(2.8) 



Here Pg is the Legendre Polynomial of the £th order and cos 7 = x • x' is the direction cosine 
between the spatial vectors defining the points x, x' . The notation r> (r<) refers to the larger 
(smaller) of r and r'. The mode function p W £ satisfies the appropriate boundary condition 
at some small value of r and q^e satisfies the appropriate boundary condition at some large 
value of r. They are both solutions to the radial equation 



ld 2 S 



h dr 2 



+ 



2 1 df 1 dh 
rh 2fh dr 2h 2 dr 



dS_ 

dr 



UJ 

7 



+ 



+ 1: 



s 



. 



(2.9) 



They satisfy the Wronskian condition 
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A WKB approximation for the modes can be derived as follows: First define 
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with 
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1 df _ f dh 
2rh dr 2rh 2 dr 

f d 2 W fld£_ f dh\ 1 dW 



hW dr 2 



3 / / 1 dW 



hdr h 2 dr J 2W dr 2h \W dr 
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u 2 + £(i + l) 
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Then solve the equation for W iteratively, the lowest order solution being W 
iteration adds two derivatives of the metric. To second order one finds 

2 



W = Q + 
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1 df f dh 



2ft \ 2rh dr 

f d 2 tt 
htt 2 dr 2 + 
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2rh 2 dr J 8fi 3 V 2rh dr 

ldf f dh\ I do 

h dr 



1 df _ f dh 
2rh 2 dr 



3f 1 (dQ\ 
h 2 drj2Q 2 dr 2hW\dr) 



(2.12) 



(2.13) 
Q. Each 
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Substitution of Eq. f[2TL]) into Eq. fl2TT0l) shows that for the WKB ansatz, C w£ = 1. 



III. COMPUTATION OF v(x,x') 



To compute an expansion for v(x,x') in powers of (x — x') one first solves Eq. ( 12.121 ) to a 
specified adiabatic order and substitutes the result into Eq. (12.81) . As shown below, the 
resulting expression can be broken into sums and integrals of the form 

roo 00 1-7/7 _|_ i \lr 

Smnp= duJCOSUj{T-r')]u 2n J2( 2i + 1 ) 

with m, n, and p integers. 



(3.1) 



Our goal is to determine the coefficient of the log(r — r') term since this is related to 
v(x,x') through Eqs. (12. 6p and (12.71) . This can be accomplished by using the Plana sum 
formula 



5>< 

n=N 



-F(N) + / dnF{n) + i 
2 Jn 



dt 



[F(N + it) - F(N - it)] 



(3.2) 



to compute the sum over i. There are superficial divergences in some of the sums over £ 
which can be isolated by expanding the summand for a given sum in inverse powers of I and 
then truncating at order l/£. If the remaining terms are subtracted from the summand then 
the result will be finite. Effectively this amounts to subtracting terms which are proportional 



7 

to 5{t — t') and its derivatives. This procedure is discussed in Ref. {^j]. After computing 
the sum over I one can expand the result in inverse powers of uj. Most of the expressions 
will be infrared divergent, however this is not important because we are only looking for the 
coefficient of the log(r — r') term. 2 Thus it suffices to put in a lower limit cutoff A on the 
integral over uj. 

The integral over uj can be computed for each term of the series. For positive powers of uj 
there are no log(r — r') terms so these do not contribute to v(x, x'). For all negative powers 
of uj there will be log(r — r') terms. For uj^ 1 one finds 

duj cos\uj{t - t')}- = -ci(A(r-r')) = -log(r-r') + ... . (3.3) 

UJ 

For all negative powers of uj, successive integrations by parts can be performed until the 
integral is in the form (13.31) . The result is 

^cosKr-rO]^ = ^ yr (r-r') 2 "log(r-r') + ... . (3.4) 

If the points are split only in the time direction then one finds that using a second order 
WKB expansion gives v to zeroth order in (t — t'), a fourth order one gives v to second order 
in (t — t'), and so forth. For Schwarzschild spacetime the lowest nonvanishing order is (t — t') 4 
which requires use of a sixth order WKB approximation. 

To compute the expansion when the points are also split in the angular direction one can 
begin by expanding the angular part of the mode functions in powers of (cos 7 — 1) so that 

P,(cos 7 ) = 1+ 1 ^ ' (cos 7 -l) + ... (3.5) 

Then the sum over i and the integral over uj can be computed as before. At the end of the 
calculation one should also expand the terms of the form (cos 7 — l)- 7 in powers of 9 — 9' and 

0-0'- 

To compute the expansion when the points are also split in the radial direction one should 
first fix r' to be either less than or greater than r. For our purpose it will not matter which 
is chosen so, as an example, we assume that r' < r. Then one expands Puii(r') in powers of 
r' — r and repeatedly uses the mode equation (12. 9p to eliminate all but first derivatives of 
Put- Thus 

Pu(r') = Pui(r) +p'ue(r)(r' - r) 



2 There is an alternative way to do the WKB expansion which gives no infrared divergences 
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+ .... 



(3.6) 



Then this expansion is substituted into the expression (12. 8ft for Ge and the WKB expansion 
is introduced to write this equation in terms of W and its derivatives as before. The terms 
will be proportional to various powers of u and £(£ + 1) multiplying either 1,1/W , or 
W'/W 2 . Those multiplying 1 are proportional to 5(r — r')5(Q — Q') and various derivatives 
of these delta functions. Therefore, they do not contribute to the calculation of v(x,x') and 
can be ignored. For the other terms one substitutes the WKB approximation for W to some 
specified order and computes the coefficients of the log(r — r') terms as before. 

To affirm the validity of the WKB scheme, we note that when this program is carried 
out for a Reissner-Nordstrom spacetime the expansion for v(x,x') agrees with the results of 
Christensen j^lj] at order (x — x') 2 . To leading order, it is also a solution to Eq. (11.21) . For 
Schwarzschild spacetime v(x, x') vanishes at order (x — x') 2 . However, by using an eighth 
order WKB approximation for the modes it has been possible to compute v(x,x') to order 
(x — x') 6 in Schwarzschild spacetime. Writing 



v(x,x')= Vij k (t - t') 2t (cos 7 - iy (r - r') 

i,j,k=0 



(3.7) 



we find 
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This expression is a solution to Eq. (11.21) to O \{x — a;') 4 ]. By that we mean that if it is 
substituted into the left hand side of Eq. (11.2 It then the resulting expression will be zero to 
0[{x-x')% 

IV. CONCLUSIONS 

We have presented a method that allows for the calculation of the part of the retarded 
Green's function v(x, x') that contributes to the tail part of the radiation reaction force for a 
massless minimally coupled scalar field in a general static spherically symmetric spacetime. 
We expect it to be straight-forward to adapt this method to the cases of the electromagnetic 
and gravitational fields in static spherically symmetric spacetimes. It may be possible to 



adapt it, for all three fields, to Kerr spacetime as we 
three fields are separable in the Kerr background 43 



1 given that the wave equations for all 



44 



45]. 



We have explicitly calculated v(x, x') in Schwarzschild spacetime to O [(x — a;') 6 ] for an 
arbitrary separation of the points. One indication of the correctness of these expressions is 
that v(x, x') satisfies Eq. ( 11.21) to the appropriate order. 

Although we have not yet computed enough terms to get an estimate of the self-force, one 



qualitative feature in our results is distinct from other related cases 51|. In the large r limit 
one can see that the leading order terms in v are all proportional to M 2 . By comparison, 
calculations by DeWitt and DeWitt and Pfenning and Poisson 23] show that when the 



spacetime curvature is everywhere small, such as is the case for a static star, then the 
leading order term is proportional to M. More specifically, in the calculation of DeWitt 
and DeWitt the metric is everywhere Schwarzshild except at the origin where they assume 
a delta function mass source. They made the approximation that the Green's function can 
be well approximated everywhere by the flat space Green's function (this is known to be 
false for regions close to the black hole event horizon) and found that the leading order term 
comes from a signal that propagates to the central condensation at r = 0, bounces off and 
comes back to the current location of the particle. 

One might be concerned that similar processes will contribute a leading order term that 
is linear in M for the black hole case. This may well turn out to be the case, and one 
should consider this factor seriously. However, before drawing any direct implications we 
caution that the case of a particle orbiting a static star is qualitatively different from that 
of a particle orbiting a black hole. For example, the calculations in Ref. 0, 2^] assume that 
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the curvature is nowhere large which is not a valid assumption for the case of a particle 
orbiting a black hole. Also, in the case of a black hole the existence of an event horizon 
precludes classical waves from scattering off the central condensation (note this is different 
from the superradiance effect occurring in the ergosphere of a Kerr black hole), even though 
such scattering can occur for orbits with r > 3M everywhere. We think it likely that this 
type of nonlocal (infrared) contribution to the self-force can be, at least approximately, 
decoupled from the quasilocal contribution we are considering here. If this is the case, then 
the nonlocal contribution can be considered separately and the result can be added to the 
quasilocal contribution to the self-force. 

The method we have developed to compute v(x,x') can be used to find the expansion 
to higher orders in (x — x'). Work is in progress to compute an expansion for v(x,x') to 
substantially higher order in Schwarzschild spacetime, which should provide a better test of 
the usefulness of the Hadamard expansion. Although the quasilocal expansion may prove 
insufficient to compute the tail part of the self-force, the facts that no regularization is 
necessary and the resulting expression is analytic make an investigation of its usefulness 
well worth the effort. Such an investigation is underway and the results will be presented 
elsewhere. 
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APPENDIX A: FIRST ERRATUM 



In the Physical Review D (PRD) version of this paper |52J there were a few notational 
and transcription errors that have been corrected in this arXiv version and also published in 
an erratum [53|. For completeness we list the corrections in that erratum here as well. 

1. In Eq. (1.1) of the PRD version the general expression for the Hadamard form of the 
retarded Green's function was given incorrectly. It should read 

G rct (x,x') = 9(x,x') !^^6[*(x,x')) - v -^le[-a{x,x')]^ (Al) 

2. The relationships between our definitions of u(x,x') and v(x,x') above and U(x,x') 
and V(x,x') of Ref. [tJ] were given incorrectly in the PRD version. To obtain the 



correct relationships first note that the convention for G ret in Ref. 52J differs from 
that of Ref. To obtain that of Ref. [7] let G re t — * Gr re t/(47r). Substituting this 
into Eq. (lAlj) and comparing with Eq. (9) of Ref. [7J gives u(x,x') = U(x,x'), and 
v (x, x') = 2V(x, x'). 

3. In Eq. (3.9) of the PRD version the expressions for the following quantities should be 
multiplied by a minus sign: u 2 oi, v in , v 103 , v 02 i, v 013 , and v 005 . 

These errors do not affect any other equations or results in the paper and the above 
corrections have been made in this arXiv version. 

We would like to thank Ardeshir Eftekharzadeh (April 2006) for first pointing out the 
errors in Eq. (3.9) of the PRD version, and Barry Wardell (April 2007) for independently 
pointing them out. P. R. A. would like to thank Hebrew University and the Department of 
Theoretical Physics at the University of Valencia for hospitality. This work was supported 
in part by the National Science Foundation under grant numbers PHY03-00710 and PHY05- 
56292. P.R.A. acknowledges the Einstein Center at Hebrew University, the Forchheimer 
Foundation, and the Spanish Ministerio de Education y Ciencia for financial support. 
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APPENDIX B: SECOND ERRATUM 



In the Physical Review D (PRD) version of this paper 52j there was a factor of 2 error 



5J. For 



that has been corrected in this arXiv version and also published in an erratum 
completeness we list the corrections in that erratum here as well. 

In Ref. 52| there was an error in the calculation of the coefficients in the Hadamard-WKB 
expansion of v (x, x'). Each coefficient in Eq. (3.9) should be multiplied by a factor or 2. A 



transcriptional error in certain coefficients was corrected in 53[]. The expansion for v is 

oo 

v(x,x') = v ijk (t - t') 2i (cos 7 - ly (r - r') k . (Bl) 

i ,j,k=0 

The corrected coefficients, including both the corrections in |53j and the factor of 2, are 



^000 = ^001 = V W0 = V 0W = V Q2 = V W l = V 011 = VQ03 = (B2) 
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